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Abstract. We discuss the transition strength between the disordered ground state 
and the basic low-lying triplet excitation for interacting dimer materials by presenting 
theoretical calculations and series expansions as well as inelastic neutron scattering 
results for the material KCUCI3. We describe in detail the features resulting from 
the presence of two differently oriented dimers per unit cell and show how energies 
and spectral weights of the resulting two modes are related to each other. We present 
results from the perturbation expansion in the interdimer interaction strength and thus 
demonstrate that the wave vector dependence of the simple dimer approximation is 
modified in higher orders. Explicit results are given in 10 th order for dimers coupled 
in ID, and in 2 nd order for dimers coupled in 3D with application to KCUCI3 and 
TlCuCl 3 . 
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1. Introduction 

Low-dimensional quantum antiferromagnets have received much interest in recent years 
since they serve as model substances allowing to investigate in detail the effects of 
quantum fluctuations and to test theoretical models. One important class of materials 
in this context consists of an assembly of dimers (two strongly coupled spins 1/2) which 
interact sufficiently weakly to guarantee that the dimer gap does not close. These 
materials are characterized by a disordered singlet ground state and a finite spin gap 
to triplet excited states. Prominent examples in this class are KCUCI3 and TICUCI3 
which have been investigated in detail in the last years by static and dynamic methods 
as well as theoretically [1]. The most detailed experimental information is obtained 
from inelastic neutron scattering (INS) experiments, which directly explore the basic 
singlet-triplet transition in all of reciprocal space [2, 3]. 

The energy of the singlet-triplet transition along the principal axis in reciprocal 
space as measured in these experiments is well described by the model of interacting 
dimers; to lowest order this is formulated as effective dimer model [4, 5], and it has been 
refined by perturbative cluster expansions up to 6 th order [3]. Here we supplement this 
analysis by discussing the dynamical structure factor. 

The dynamical structure factor for spins localized on a Bravais lattice is defined as 

/oo 
dt e - <wt (5 a (q,t)5' J (-q,0)> (1) 
-00 

where 

^(q,t)=]Te- iqR ^(R,t) (2) 

R 

is the Fourier transformation of the spin operators at lattice sites R. The superscripts 
a, (3 denote the spin components and the brakets (• • •) thermal expectation values (which 
for T = reduce to groundstate expectation values (0|---|0)). Apart from known 
prefactors, Eq. (1) reflects the spectral weight from the magnetic neutron scattering 
cross section [6]. 

If we consider transitions from the ground state |0) to some well-defined eigenstate 
\n) with energy cu n (q), we obtain 5-peaked contributions to the dynamical structure 
factor 

S^(q,a;) = ^(0|^(q)|n)(n|^(-q)|0) S(u - u B (q)) (3) 

n 

= £/^(q)%-c„(q)). (4) 

n 

In interacting dimer materials, INS probes directly the transition from the (singlet) 
ground state |0) to the lowest (triplet) excitation \t) and we will reduce our discussion 
to this contribution to the dynamical structure factor. Owing to the rotational symmetry 
of the underlying Heisenberg model it is sufficient to calculate Jf z (q) only and we use the 
shorthand I sm (q) := i^*(q) to denote the lowest triplet (single magnon) contribution to 
the spectral weight. 
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The INS investigation of the materials KCUCI3 and TICUCI3 at finite magnetic 
fields provides direct verification of this point as reported in [7]. 

The discussion of our results is organized as follows. In section II we will give 
theoretical results for a ID array of interacting dimers. This model was treated before [8], 
it is, however, instructive to demonstrate for the simple ID case, that existing standard 
expansions are modified by additional terms which emerge starting in second order. In 
addition we present the quantitative changes for the transition strength comparing first 
order to 10 th order results to show the effect of high order calculations. In section III we 
discuss interacting dimers in a 3D network by presenting in parallel neutron scattering 
results for the material KCUCI3 and series expansions to 2 nd order. The same type of 
additional terms as in ID is obtained in this calculation and corrects the results for 
the dynamical structure factor as obtained in the random phase approximation (RPA) 
before, see Ref. [9]. These RPA results are found to be correct only to first order. 
Section IV gives our conclusions. 

2. Alternating chain 

First we consider the one dimensional (ID) alternating S — 1/2 spin chain with isotropic 
nearest-neigbour interactions. The hamiltonian of this model is of the following form 

N 

H = J^(S 1 (n)-S 2 (n) + AS 2 (n)-S 1 (n + l)), J > 0. (5) 

n=l 

Here, the alternating chain is described as a system with N unit cells with two spins 
each and periodic boundary conditions are used. There are two exchange constants, J 
and AJ, for A = the ground state of the system consists of singlets on the intracell 
bonds (n, 1) — (n, 2). These local singlets can be excited to triplets which remain gapped 
excitations when switching on A, < A < 1. In the limit A = 1 we arrive at the well 
known Heisenberg antiferromagnet (HAFM) with pairs of S = 1/2 spinons as lowest 
gapless excitations. Other related models are described in Ref. [10]. 

The triplet excitation energies u(q) have been obtained by perturbation expansion 
in A up to 9 th order in Ref. [8] and to 10 th order in Ref. [5] using the cluster expansion 
approach. 

2.1. The dynamical structure factor 

Turning to the calculation of the structure factor for a system with hamiltonian (5), 
we note (see Eqs. (1,2)) that for this calculation we have to specify the positions of the 
spins, R, in space (whereas the eigenvalues depend solely on the exchange constants). 
In slight generalization of a strictly linear geometry we allow for our calculation the 
separation d of two spins in one unit cell to be different in magnitude and direction 
from the separation a of two adjacent spins in different unit cells (note that a defines 
the overall chain direction). The resulting geometry is shown in Figure Al and makes 
clear the relation to the real 3D systems to be dealt with in the next section: the 
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ID chain denned in Eq. (5) can alternatively be looked at as a ladder with rung and 
diagonal interactions only. For the chain geometry shown in Figure Al the spectral 
intensity up to first order in A is obtained as follows 

I sm (q) = sin 2 ^(l + ^Acos(qa)) + 0(A 2 ). (6) 

Here, q is the wave vector, d and a the separation of the spin sites within and between 
the dimers, respectively We note that the term oc sin 2 (qd/2) is typical for systems 
consisting of isolated dimers, it is known as the dimer structure factor [11]. The first 
order correction in A adds an additional modulation to the intensity, which depends on 
the ratio a = ||d||/||a||. 

Using the cluster expansion method (see Appendix A) we have systematically 
calculated the series in A for the intensity up to the tenth order. This requires linked 
clusters consisting of maximum 10 bonds. The resulting series can be split into three 
different terms: 

7 sm (q) = 5 c (q, A) + B s (q, A) + A(q, A). (7) 
To illustrate the result we give the series up to fourth order!:: 

B c (q, A) = sin 2 ^- ^ ^ cos(jqa), 

3=0 
4 

B s (q, A) = sin qd ^ Vj sin(jqa) 



(8) 



where 

,o = 1-^A 2 -Aa3 + i |^, , o = 0, (9) 

a - \ V 5 A 3 + 41 A 4 v + 7 A 3 131 X 4 

^-2 A "8 A "192 A + 2304 A ' ^ ~ 8~ A + 192~ A " 4608 A ' 

^ 2 ~ 16 A + 48 A + 1024 A ' U2 ~ 96 A + 4608 A ' 

5 x 3 155 4 23 4 

^ ~ 64 A + 2304 A ' U3 ~ 2304 A ' 

35 x 4 n 
^= 1024 A ' Ui = °' 

and 

A(q,A) = T ^A 4 (cos(2qd) -cos(2qa)). (10) 

The terms in B c (q, A) are consistent with previous publications [8], whereas B s (q,X) 
and A(q, A) contain additional corrections. They originate from a complete expansion 
of both the ground state |0) and the first excited state \t). If one assumes that 
(0\S2i\t) = — (0\S2i+i\t) for the matrix elements on even and odd sites, one ends up 
with B c (q, A) only. However this is only correct to first order and in general we have 

| Higher order terms are available on request. 




Figure 1. Structure factor with A = 0.4, 0.8, a || d and ratio d/a = 10/20. 



(0\S2i\t) 7^ —(0\S2i+i\t). This inequality arises from virtual states with odd parity under 
exchange of two triplets which occur during the perturbation expansion for the first time 
in second order. 

In Figure 1 and Figure 2 we show some typical plots of the intensity I sm (q) for two 
different ratios a = 10/15, 10/20 and two different coupling strength A = 0.4, 0.8 and 
strictly linear geometry, d||a (then, only the component of wavevector in chain direction 
enters). The difference between the zeroth and the first order emerge very clearly. The 
higher order terms emphasize the modulation originating from the two length scales ||a|| 
and lldll. 



2.2. Sum rule 

The total integrated scattering intensity has a well defined magnitude, determined by 
the local spin length through the following sum rule: 

For the one dimensional alternating chain the contribution from the one-magnon part to 
the total spectral weight is calculated from eq. (7); the integral reduces to the constant 
part fiQ of (9), since only non oscillating terms survive the integration: 

J sm = -fl--A 2 --A 3 + ^A 4 + ...^) <- S=-, A « 1.(12) 
4 I 16 32 1536 ~ A 2' v ; 




Figure 2. Structure factor with A = 0.4, 0.8, a || d and ratio d/a = 10/15. 



For the noninteracting case (A = 0) the sum rule is exhausted by the one triplet 
excitation since it is an exact eigenstate: From (3) we see that this excitation gives 
the only no n- vanishing matrix element. Switching on the coupling between the dimers, 
more and more intensity goes in two or more magnon scattering processes. A theoretical 
discussion of the multimagnon states for the one dimensional alternating chain is given 
in Refs. [8, 12, 13]. 



3. The 3D dimer substances KCuCl 3 and TlCuCl 3 



In this section we extend the calculation of singlet-triplet intensities to three dimensional 
substances such as KCuCl 3 and TlCuCl 3 . These compounds are weakly interacting 
quantum spin systems which exhibit an excitation gap. Similar to the alternating chain 
discussed in the previous section this is based on the existence of strongly interacting 
bonds forming dimers. In these materials the orientation af the dimers alternates, i. e. 
each unit cell consists of four spins with two differently oriented dimers. 

The following considerations are valid for dimer systems composed of two dimers 
per unit cell. We write for the Fourier transformed spin operators 

2 



S 2 (q) = 



1 



f2N 



EE 



-jq(n+R fc ) 



k=l 



^> + R fe + |) + e^n + R fc -|; 



(13) 



where n denotes the unit cell, the center of the dimer and d^ the separation of the 
two spins forming the dimer. The first sum in Eq. (13) is taken over all unit cells. In the 
case of weakly interacting dimers the localized triplet states are replaced by Bloch-like 
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Table 1. Considered interactions in KCUCI3 (and TICUCI3). Vectors gk denote the 
distances between the dimer centers. 



distance gk Interactions between 

between dimers equivalent spins (i = j) non equivalent spins (i 7^ j) 



gi = a 




^(100) 


J (100) 


g2 = 2a 


f c 




Jl 

J (201) 


g3 = a + 


i(b + c) 




Jl 


g4 = a - 


i(b-c) 




ji 



Table 2. Values of exchange interactions in KCuCi3[5] and TlCuCl3[3]. The 
intradimer interaction J is given in meV, intcrdimer interactions are given in units 



of J for the respecti 


ve compounc 


1 


exchange constant 


KCuCl 3 


TlCuCl 3 


J 


4.25 meV 


5.68 meV 


^(100) 


0.00 


0.06 


j' 


0.10 


0.30 


ji 

J (201) 


0.18 


0.45 




0.20 


0.16 


Jl 


0.05 


-0.10 



triplet modes which propagate due to the interaction network between the dimer units. 
The details of the explicit calculation of the transition matrix elements for a Bravais 
dimer lattice are illustrated in Ref. [8]. 

In Tab. 1 we describe the interaction network listing the basic lattice vectors 
associated with nonzero interdimer interactions for the materials considered. In Tab. 2 
we give the numerical values of the intradimer exchange J (in meV) and of the interdimer 
exchange interactions (in relative units) for the compounds KCuCl 3 and TlCuCl 3 . 
(Slightly improved values for KCUCI3 have been determined in [17], the difference, 
however, is not visible in the Figures 3 and 5 below.) 

The different dimer orientations in the unit cell do not affect the dimer lattice 
directly, i. e. the lowest excitation does not depend on the dimer orientation. But the 
full translational symmetry is obviously reduced if the dimer sites are distinct by their 
orientation. 

In analogy to phonons in a lattice with a basis there will be two excitation modes. 
Therefore we use the following zeroth order ansatz for the one triplet wave function, 
which manifest translational symmetry: 

|q) (0) = ~J=Y, E Cfce " q(n+Rfc) l n + ^) {0) - (14) 

The states |n + R fc )(°) denote a triplet at site n + R fe with all the other sites occupied by 
singlets. We have introduced the coefficient c k to take into account the different dimer 
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orientations. The c^'s are determined requiring that ci|n + Ri) + C2|n + R 2 ) is diagonal 
in the subspace of the one triplet excitations. Two solutions (±) are obtained: 

(c+ 4) = (1, 1) or (cr, c 2 -) = (1, -1). (15) 

The solution c~ is connected to umklapp scattering where q — > q + u. Umklapp 
processes are possible in KCuCl 3 and TlCuCl 3 with an integer number of the reciprocal 
lattice vector u = b* or u = c*. 

The resulting energies are denoted by w ± (q) where we describe both modes in 
the first cristallographic Brillouin zone. Up to first order we obtain the well known 
dispersion relation [5, 4] in units of J: 

2 4 

2^Acos( gi q) + 0(A 2 ). (16) 

i=l i=3 

Here and in the following we use some short-hands of various combinations of coupling 
constants (e^ and 7^ will be needed for the intensity calculation below): 





= |(2J(ioo) - 


■^(100))) 


^1 = 


-(2J(ioo) + ^(100) )> 


7i = 


4 J (ioo)> 


@2 


- 4 J (201)> 




£2 = 


4 J (201)> 


72 = 


4 J (201)' 


ft 




j ;iii))' 


£3 = 


4(^(111) + J ('lII)). 


7? = 


4/(±^li§) + ^lII) 


Pa 






e 4 = 


£3, 


7f = 


73 ± 



To obtain the transition matrix element one has to expand both the ground state and 
the one triplet state perturbatively in the coupling constants. To indicate the order we 
take all interdimer couplings proportional to a constant A. 



3.1. The ground state 

There are four different directions gj in which we find interdimer interactions (see 
Tab. 1). As well as in the one dimensional case the unperturbed ground state is a 
product of singlets placed on the rungs: 

|G> ( ° ) =ni a n+Ri>l«n + H a > = |5>. (18) 
n 

Due to the different orientations we distinguish between singlets at n + Ri and n + R 2 . 
Up to the second order the ground state including all relevant states for the structure 
factor is 

|G)( 2 ) = a \S) + ^ J2 E E & (1 + ^)l n + R *' n+R fe+& >(°< ) (19) 

i=l k=l n 
/o 4 2 

- \ E EE^i n+R ^ n + R ^+&+^) (0,0) ( 2 °) 

i,j=l k=l n 

- ^ ttE^>+R tl n+R,+g i -g J ) (0 ' 0) (21) 

»,j=l k=l n 
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+ states with three or four triplet excitations. (22) 

The indices i,j are linked to the different interaction directions and k counts the two 
dimer sites in the elementary cell. We further denote the states having two triplets at 
sites r and r' with well-defined S tot and S£ ot by |r, r')( Stot ' s *<*\ a is a normalization 
factor which guarantees that (G\G) {2) = 1 + C(A 3 ): 

a = l- 3 -N^(3l (23) 
i=i 

Note that N labels the number of unit cells, whereas 2N is the number of dimers in the 
system. 



3.2. One triplet excitation 

The expansion of the wave function for the one triplet excitation in the interdimer 
interactions is obtained to first order by acting with the Hamiltonian Hi on the state (14) 
and some subsequent normalization. In addition to simple propagation of the triplets 
this leads to the generation of two triplet excitations |r, r')( 1,0 ) where (<St t, Stot) = (1> 0) 
is the total spin and total magnetization and r, r' label the lattice sites occupied by 
triplets. There are also three triplet excitations with the same spin quantum numbers. 
Second order terms contribute to the third order§ of the intensity only and will not 
be calculated here. However, normalization of the wave function has to be done up to 
second order terms. 



3.3. The dynamical structure factor 

To leading order one expects a dimer-like structure factor as in Sec. 2. In fact there are 
two different contribution due to the two dimer structure in the elementary cell: 

2 



/±(q) = £>±(q) + ^(A 1 ) 



. qdi qd 2 

sm ± sm 

2 2 



+ 0(\ 1 ). (24) 



The indices ± refer to (16) and correspond to symmetric resp. antisymmetric modes for 
q in the first Brillouin zone. However, the role of symmetric and antisymmetric modes 
is interchanged for q — > q + r where e "~( R 2-Ri) — _i From now on q stays in the first 
Brillouin zone. 

The representation of (24) is instructive in order to emphasize that in general there 
are contributions from two different modes for a wavevector q in the crystallographic 
Brillouin zone. The total contribution 

I + (q) + /_ (q) = 2 sin 2 ^ + 2 sin 2 ^ (25) 

reproduces correctly the structure factor J du>S(q,u) as calculated from (1) in lowest 
order. 

§ All contributing terms can be obtained on request. 
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A finite contribution for the excitation mode with energy u (q) requires qd x 7^ qd 2 . 
Taking into account the dimer orientations d 12 where 

di = 0.48a + 0.10b + 0.32c, 

d 2 = 0.48a- 0.10b + 0.32c (26) 

we deduce the condition qt 7^ 0, in agreement with the experimental observations. 

Considering higher order corrections to the ground state and the first excited state 
as presented below, we get the following result valid up to second order: 

7±(q) = |^±(q)(l " ^±(0) + ^ (q) - S ± (q) - S^q)) 

+ iD ± (q)(cos^i±cos^)A ± (q) (27) 

where 

2 4 

fi±(q) = 2^Acos(qgi) ±2^Acos(q gi ), ( 28 ) 

i=l i=3 

4 2 4 

S±(q) = ^ [3/3, 2 - 01 cos(2q gl )] + 2 ^ /3 i£i cos(q & ) ± 2 ^ # £i cos(q & ), (29) 

i=l i=l i=3 

2 4 

S x (q) =4^7- [l + cos(qg i )] + 2^ [ 7 f + 7"' + ^lt% cos(q gl )] , (30) 

i=l i=3 
2 4 

A±(q) = 2^A 7j sin(q & ) + 2^A( 7 r ±7 l + )sin(qg J ) (31) 

i=l i=3 

The dynamical structure factor contains contributions from both excitation modes: 

5"(q, w) = 7 + (q) <5(u; - W +(q)) + 7_(q) *(a; - uT(q)) , (32) 

where w ± (q) denotes the one triplet excitation energy as in (16). We note that fi(q) is 
the energy of the one triplet excitation w + (q) up to first order. If we neglect the terms 
Sj(q) and A±(q) the result reduces to RPA-like calculations [9] where 7±(q) oc l/o; ± (q). 

In order to demonstrate the effect of interdimer interactions on the dynamical 
structure factor we show in Figure 3 theoretical results for KCuCl 3 (top panel) and 
TlCuCl 3 (bottom panel) in th (noninteracting dimers), 1 st and 2 nd order, exchange 
parameters are taken from Tab. 2. As in the INS experiments to be discussed in the 
next subsection and as in Ref. [9], the variation of the spectral weight of the triplet 
excitation I + + 7_ with wave vector is shown along the (0,x,x) direction of reciprocal 
space, such that both modes contribute with finite weight. Evidently, higher order 
corrections are more important for TICUCI3 with its larger exchange constants, but even 
for these larger values the comparison of different orders seems to indicate convergence. 
Exhaustive experimental results for both 7 + and 7_ are available for KCUCI3 and will be 
discussed in the next subsection. Clearly, devoted INS experiments for TlCuCl 3 are of 
considerable interest (see also [4] and [3]); the ratio of spectral weights at wave vectors 
corresponding to maximum and minimum intensity appears to be a reasonable quantity 
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Figure 3. Total spectral weight i+ + 7_ along the (0, a;, a;) direction in th , 1 st and 
2 nd QldeT Top panel . KCuCl3j bottom p ane i : TlCuCl 3 . 



to test the agreement with our theoretical results. Present results do not indicate that 
higher than 2 nd order calculations are required. 

3.4- Experiment 

The INS results on the material KCuCl 3 were collected at the IN3 neutron spectrometer, 
Institut Laue-Langevin, Grenoble (France). Standard focussing geometry was adopted 
for all energy scans performed under constant final energy Ef = 13.7meV. A pyrolitic 
graphite (PG) filter in front of the analyser was further used to suppress higher order 
contaminations. 
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Figure 4. Typical fits of the neutron profiles for wave vectors q taken at selected 
(0, x, x) values in (r.l.u) 



The INS profiles were obtained at fixed T = 2K for the wave vectors along 
the (0, x, x) direction of reciprocal space, which is suited to demonstrate the issues 
introduced above. For this purpose, a KCuCl 3 single crystal was aligned for scattering 
in the b*c* plane. The spectral weight of the triplet excitation was determined from 
global least squares fits to the measured neutron profiles, assuming gaussian peaks on 
top of a common background. The center of the peaks was further fixed at the energies 
cu ± (q) resulting from the analysis presented in Ref. [5]. (see Tab. 2) Our present results 
complete the experimental investigation of the b*c* plane summarized in Ref. [4], and 
references therein. 

In Figure 4 typical fits of the neutron profiles are shown for wave vectors q at 
selected (0,x, x) values, in reciprocal lattice units (r.l.u.) of the unit cell. In accordance 
with the theoretical expectations, both excitation modes o; ± (q) are visible along (0, x, x) 
but the spectral weight l±(q) strongly depends on x. Continuous lines denote the 
global least squares fit function, symbols the profiles in neutron counts. The statistical 
tolerance scales according to the neutron counts. In Figure 5, the fitted spectral weight is 
compared to the model expectations previously introduced. In the top panel, the total 
spectral weight F 2 (I + + J_) (full line) is compared to the experimental observations 
(symbols). The only free parameter is an overall scaling factor accounting for the size 
of the sample, both the plain calculation (dashed dotted line, second order) and the 
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wave vector (0, x, x) [r.l.u] 




wave vector (0, x, x) [r.l.u] 



Figure 5. Spectral weight of the triplet excitation for KCuCbj along the (0,x,x) 
direction in reciprocal space. Symbols indicate experimental data from least squares 
fit to the profiles. Top panel: theoretical result for F 2 {I + + 7_) in second order (full 
line) and for noninteracting dimers (dashed line), F 2 (dotted line) and (7+ + 7_) in 
second order (dash-dotted line). Bottom panel, full line: theoretical result in second 
order for the relative spectral weight. 

squared magnetic form factor F 2 (dotted line) [14] are shown separately for convenience. 
In the bottom panel, the relative spectral weight + 7\_)(q) is compared to the 

experimental observations, as indicated. The graphical representation underlines the 
redistribution of the spectral weight among I + and J_ which occurs along the (0, x, x) 
direction of reciprocal space. 

From Figure 5, reasonable agreement between predictions from the dimer model and 
experimental results is concluded. The spectral weight is dominated by the bipartite 
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dimer structure, which governs the result in the noninteracting dimer limit, but the 
existence of higher order corrections is clearly seen close to the second maximum (x ~ 3). 
Details of higher order corrections remain almost beyond statistics for KCUCI3 but 
may become more pronounced in the sister material TICUCI3. Our results improve 
on the previous RPA calculations [9] which are correct only to first order. Related 
investigations along different directions of the b*c* plane were successfully compared to 
RPA calculations in Ref. [15]. The relative spectral weight (bottom panel of Figure 5) is 
very well described already in the noninteracting dimer picture (not shown in Figure 5), 
higher order corrections are below statistical significance for this quantity. 



3.5. Sum rule 



We calculate the one magnon contribution to the total integrated scattering intensity in 
order to check the sum rule (12). As seen in (32) the dynamical structure factor consists 
of two parts. Integrating over q, only non oscillating terms survive, giving 



T — 



4 



1 " E # - 4(7 2 + ll) - 2(1 f + if + if + 11 



+C(A 3 ).(33) 



Using the coupling constants as calculated in Ref. [5] we estimate here the intensity 
which goes into higher order scattering processes like two or more magnon scattering. 
We obtain T sm = 0.7217 which means that 96.23% of the total scattering intensity is 
concentrated in the lowest triplet excitation. Although the interactions in TlCuCl 3 are 
more pronounced most of the scattering intensity still goes in the one magnon process 
which is reflected by inserting calculated coupling constants [3] in (33): T sm = 0.7021 
or 93.62% respectively. The absolute experimental determination of the spectral weight 
from dimers is exemplified in Ref. [16], but a devoted investigation of the materials 
KCUCI3 and TICUCI3 has not been performed up to now. 



4. Conclusion 



We have presented series expansions for the dynamical structure factor valid generally 
for lattices with two dimers per unit cell. Applying our results to the interacting dimer 
material KCUCI3, we have shown that results obtained by inelastic neutron scattering 
are reasonably well described by the theoretical calculations. Our expressions apply 
as well to the sister material TICUCI3, which shares with KCUCI3 the structure of the 
exchange couplings, but has larger exchange strengths. For the specific direction in 
q-space considered here our results show that higher order terms are not relevant for 
relative spectral weights (see Figure 5 as measured in KCUCI3 with present intensity) 
and we expect that this is generally true. Second order shifts, however, show up in 
absolute spectral weights [17], most clearly in TlCuCl 3 . 

The materials KCuCl 3 and TlCuCl 3 recently have been demonstrated to undergo 
field-induced magnetic ordering. The evolution of the excitation modes at finite 
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Figure Al. Linear arrangement of dimers indicating the similarity of alternating 
chain (interdimer interactions only between ^(n) and S\(n+ 1)) and ladder (general 
interdimcr interactions) geometries 



magnetic field has been described in a comprehensive theoretical study [18], albeit 
limited to the energy of the excitations. Theoretical investigations of the spectral weight 
in an external magnetic field are now under preparation. 

Appendix A. Cluster Expansion 

In this appendix we briefly summarize the method of cluster expansion for the dynamical 
structure factor in the case of the alternating chain. Some detailed considerations can 
be found elsewhere [19]. As shown in Figure Al the cristallographic unit cell contains 
two spins. Thus the Fourier transform of the spin operator splits into two parts and 
reads as: 

S a (d) = e ~ iqna (e - *^^ (n) + e^S«{n)} . (A.l) 

n 

As before ||a|| = a is the distance between neighbouring spins and d denotes the spacing 
between the two spins on a dimer. In our notation q is the projection of the wave vector 
q on the chain direction. 

Using translational invariance with respect to the center of the dimer we obtain for 
the singlet-triplet transition amplitude: 

Wq) = E ^ W + + e^SM + e " qd ^ (n)} (A.2) 

n 

where A$(n) = (0\S*(0)\t)(t\S?(n)\0), i,j = 1,2. (A.3) 

Here, the sum is taken over all integer numbers n. However, it is more convenient to 
calculate the functions A\Hn) for positive numbers n. This is feasible making use of 
inversion symmetry wrt to the dimer center, implying 

AlK-n) = A%(n) and Ag(-n) = A%(n). (A.4) 

Inserting (A.4) into (A.2) one arrives at the following result: 

J sm (q) = 2 J2 ( A Ti( n ) + A 22 W) cos(qna) (A.5) 

n>0 

+ 2 E ( A i2 ( n ) cos(qna - qd) + A^(n) cos(qna + qd)) (A.6) 

n>0 
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+ A™(0) + A%(0) + e iqd A^(0) + e- iqd A^(0). (A.7) 

Now, functions A^(n) enter for positive n only. In the limit of noninteracting dimers 
only the terms with n = (last line) survive. 

At first glance the functions Afj(n) are groundstate expectation values which can 
be computed by the well established cluster expansion method [20]. The projection 
operator V = \t)(t\ has to be evaluated from the one magnon states where i labels 

the lattice site. By means of degenerate cluster expansion these states are generated 
order by order [21]. Then we find for projection operator: 

^ = |i><i| = E^" 1 )^ (0 >^ (, " ) i- (A- 8 ) 

ij 

g is the overlapping matrix of the 

9a = (^ w ). (A.9) 

To invert g we use the fact that g is the unit matrix for A — > 0: 

9 = 1 + 9- (A.10) 
Owing to the matrix norm \\g\\ < 1 we apply a geometric series to invert g: 

oo 

{I+~g)- l = Y.~9 l - ( A - n ) 

i=0 

Now we have everything at hand to calculate the singlet-triplet-intensity of the 
dynamical structure factor: Apply degenerate perturbation theory to obtain the states 
and V . Then calculate g and invert this matrix by using (A. 11). Finally, apply 
non degenerate perturbation theory to compute the functions A\*{n). 
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